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Matrix ansatz and combinatorics of the /c-species PASEP 

Olya Mandelshtam 

Abstract 

We study a generalization of the partially asymmetric exclusion process (PASEP) in which 
there are k species of particles of varying weights hopping right and left on a one-dimensional 
lattice of n sites with open boundaries. In this process, only the heaviest particle type can enter 
on the left of the lattice and exit from the right of the lattice. In the bulk, two adjacent particles 
of different weights can swap places. We prove a Matrix Ansatz for this model, in which different 
rates for the swaps are allowed. Based on this Matrix Ansatz, we define a combinatorial object 
which we call a fc-rhombic alternative tableau, which we use to give formulas for the steady 
state probabilities of the states of this fc-species PASEP. We also describe a Markov chain on 
the 2-rhombic alternative tableaux that projects to the 2-species PASEP. 
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1 Introduction 

The (single-species) partially asymmetric exclusion process (PASEP) is an important non-equilibrium 
model that has generated interest in many areas of mathematics. Partly this is due to the existence 
of an exact solution (i.e. there are explicit formulas) for the stationary distribution for this process, 
which makes it a useful example in the study of non-equilibrium processes. Moreover, the PASEP 
has a rich combinatorial structure. Most notably, Corteel and Williams gave a beautiful com¬ 
binatorial interpretation of the stationary distribution of the PASEP in terms of certain tableaux 
called permutation tableaux, which are in bijection with several related objects, such as alternative 
tableaux, staircase tableaux, tree-like tableaux. 

The single-species PASEP describes the dynamics of particles hopping on a finite one-dimensional 
lattice on n sites with open boundaries, with the rule that there is at most one particle in a site, 
and at most one particle hops at a time. Figure shows the parameters of this process, with the 
Greek letters denoting the rates of the hopping particles. 

More precisely, if we denote the particle by d and the hole (or absence of a particle) by e, and let X 
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Figure 1: The parameters of the PASEP. 


and V be any words in {e, d} then the transitions of this process are: 

XdeY ^ XedY eX ^ dX Xd ^ Xe 

Q 


where hy X ^Y we mean that the transition from X to Y has probability n being the length 
of X (and also Y). 


The fe-species PASEP is a generalization of the single-species PASEP, where we now have k particle 
types of varying weight hopping on a finite one-dimensional lattice of n sites. We call the heaviest 
particle type d, the hole (or absence of any particles) type e, and the rest of the particle types 
from heaviest to lightest, type ai through type ak-i- A heavier particle type can swap places with 
a lighter particle type with rate 1 if the heavier particle is on the left, and rate qi if the heavier 
particle is on the right, where Qi is fixed to be some parameter out of a list of parameters {(?*} (see 
Section]^ for the precise rules). Like in the single-species PASEP, particles of type d can enter on 
the left of the lattice and exit on the right of the lattice. Note that for qi = q for all i, when k = 1, 
we recover the original single-species PASEP, and when k = 2, we obtain the two-species PASEP, 
which has been studied, among others, by Uchiyama |13| . Schaeffer [^, and Ayyer |1]. 


Derrida et. al. |5] provided a Matrix Ansatz solution for the stationary distribution of the single¬ 
species PASEP. The Matrix Ansatz is a theorem that expresses the steady state probabilities of a 
process in terms of a certain matrix product for some matrices that satisfy some conditions that 
are determined by the process. Uchiyama provided a similar Matrix Ansatz solution, along with 
corresponding matrices, for the two-species PASEP. In this paper, we extend the result of Uchiyama 


to a Matrix Ansatz for a type of A:-species PASEP in Theorem 5.1 





Figure 2: (a) an alternative tableau, (b) a RAT, and (c) a 3-RAT 

The tableaux of Corteel and Williams that provided a combinatorial interpretation for the stationary 
distribution of the single-species PASEP were generalized in |9] to certain tableaux called rhombic 
alternative tableaux (RAT) for the two-species PASEP. In this paper, we present even more general 
tableaux called the fc-rhombic alternative tableaux (/c-RAT), that give a combinatorial interpretation 
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in Theorem 6.7 for the stationary distribution of the /c-species PASEP. Figure shows examples of 
the alternative tableaux, the rhombic alternative tableaux, and the /c-rhombic alternative tableaux. 
Note that 2-RAT are the same as RAT, and 1-RAT are the original alternative tableaux. 


A common strategy to prove the combinatorial formula for the tableaux is to use the Matrix Ansatz, 
as in |3|. Another interesting and more direct proof of this formula is to explicitly construct a 
Markov chain on the tableaux, and show that it projects to the Markov chain of the underlying 
process, as in |2]- As our final result, we construct such a Markov chain on the RAT that projects 
to the two-species PASEP. 

Our paper is organized as follows. In Section we describe the two-species PASEP. In Section of 
this paper, we provide a proof for Theorem |2.9| by explicitly defining the matrices that both provide 
the weight generating function of the RAT, and also satisfy the Matrix Ansatz hypothesis. 


In the second half of the paper, we describe a generalization of the two-species PASEP process to 
a /c-species PASEP, and a generalization of the RAT to the A:-RAT. Our proofs are analogous to 
the proofs for the two-species case. In Section of this paper, we introduce the A:-species PASEP. 
In Section we provide a Matrix Ansatz theorem that describes the stationary probabilities of 
this process as certain matrix products. In Section we define the /c-RAT, which provides an 
interpretation for the stationary probabilities of the fe-species PASEP. Finally, in Section we 
describe a Markov chain on the RAT that projects to the PASEP. 


Acknowledgement. The author gratefully acknowledges Lauren Williams for her mentorship, and 
also Sylvie Corteel and Xavier Viennot for many fruitful conversations. The author also thanks the 
France-Berkeley Fund and the NSF grant DMS-1049513 that supported this work. 


2 Previous results on the two-species PASEP 


First we describe the two-species PASEP and the associated rhombic alternative tableaux to motivate 
the more general /s-species process that we study in this paper. 


The two-species partially asymmetric exclusion process (PASEP) has been studied extensively as an 
interesting generalization to the single-species PASEP. The two-species PASEP has two species of 
particles, one “heavy” and one “light”. The “heavy” particle can enter the lattice on the left with rate 
a, and exit the lattice on the right with rate (3. Moreover, the “heavy” particle can swap places with 
both the hole and the “light” particle when they are adjacent, and the “light” particle can swap places 
with the hole. Each of these possible swaps occur at rate 1 when the heavier particle is to the left of 
the lighter one, and at rate q when the heavier particle is to the right (we simplify our notation by 
treating the hole as a third type of “particle”). The parameters of the two-species PASEP are shown 
in Figure where the • represents the “heavy” particle, and the • represents the “light” particle. 
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Figure 3: The parameters of the two-species PASEP. 

More precisely, if we denote the “heavy” particle by d, the “light” particle by a, and the hole by e, 
and let X and Y be any words in {d, a, e} then the transitions of this process are: 

XdeY ^ XedY XaeY ^ XeaY XdaY ^ XadY 
Q q Q 
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eX dX 


Xd^Xe 


where hy X we mean that the transition from X to y has probability n being the length 
of X (and also Y). 

Notice that since only the “heavy” particle can enter or exit the lattice, the number of “light” particles 
must stay fixed. In particular, if we fix the number of “light” particles to be 0, we recover the original 
PASEP. 

Uchiyama provides a Matrix Ansatz along with matrices that satisfy the conditions, to express the 
stationary probabilities of the 2-species PASEP as a certain matrix product. 

Theorem 2.1 ([l3]). Let W = Wi ... Wn with Wi G {d, a, e} for 1 < i < n represent a state of the 
two-species PASEP of length n with r “light” particles. Suppose there are matrices D, E, and A and 
vectors (tcl and |u) which satisfy the following conditions 


DE = D + E + qED DA = A + qAD AE = A + qEA 


{w\E = —{w\ 


a 


Dlv) = -\v) 


then 


1 " 

Prob(iy) = ——(rcl '^(Wi=d) +^l(Wi=a) +El(^Wi=e) 1^^) 

^ti.r 


2=1 


where Zn,r is the coefficient of y” in (^l(-P+K^+-^) h) _ 


{w\A'^\v) 

This result generalizes a previous Matrix Ansatz solution for the regular PASEP of Derrida et. al. 

in |5]. 


In a previous paper, Mandelshtam and Viennot j9] introduced the rhombic alternative tableaux 
(RAT) which generalize the alternative tableaux and provide an explicit combinatorial formula for 
the stationary probabilities of the two-species PASEP. We describe the rhombic alternative tableaux 
below. 




Figure 4: The rhombic diagram Figure 5: An example of a RAT of size (9,3,4) 

T{daaddedaee). with type daaddedae and weight a^jd^q^. 


Definition 2.2. Let LP be a word in the letters {d,a,e} with k d’s, i e’s, and r a’s of total length 
n := k -I- £ -h r. Define Pi to be the path obtained by reading W from left to right and drawing a 
south edge for a d, a west edge for an e, and a southwest edge for an a. (From here on, we call any 
south edge a d-edge, any east edge an e-edge, and any southwest edge an a-edge.) Define P 2 to be 


4 






























the path obtained by drawing i west edges followed by r southwest edges, followed by k south edges. 
A rhombic diagram r(Vl^) of type W is a closed shape on the triangular lattice that is identified 
with the region obtained by joining the northeast and southwest endpoints of the paths Pi and P 2 
(see Figure]^. 



Figure 6: The tiles de, da, and ae. Figure 7: (a) maximal and (b) minimal hexagons. 

Definition 2.3. A tiling P of a rhombic diagram is a collection of open regions of the following 
three parallelogram shapes as seen in Figure the closure of which covers the diagram: 

• A parallelogram with south and west edges which we call a de tile. 

• A parallelogram with southwest and west edges which we call an ae tile. 

• A parallelogram with south and southwest edges which we call a da tile. 

We define the area of a tiling to be the total number of tiles it contains. 

Definition 2.4. The size of a RAT of type W is (n,k,r), where k is the number of d's in W, r is 
the number of a’s in W, and n is the total number of letters in W. 

Definition 2.5. A e-strip on a rhombic diagram with a tiling is a maximal strip of adjacent tiles 
of types de or ae, where the edge of adjacency is always an e-edge. A d-strip is a maximal strip 
of adjacent tiles of types de or da, where the edge of adjacency is always a d-edge. An a-strip is a 
maximal strip of adjacent tiles of types da or ae, where the edge of adjacency is always an a-edge. 

Definition 2.6. To compute the weight wt(P) of a filling F, first a q is placed in every empty tile 
that does not have an a below it in the same e-strip or a /3 to its right in the same d-strip. Next, 
wt(P) is the product of all the symbols inside F times for F a filling of size {k + i + r,r, k). 

An example of a RAT is shown on the left of Figure]^ 

It is easy to show with a weight-preserving bijection that we define in terms of certain “flips” on the 
tilings, that the sum of the weights of the tilings of r(VF) is independent of the tilingQ 

Proposition 2.7 ([^ Proposition 2.8). Let W be a word in {d, a,e}. Let 7i and F 2 represent two 
different tilings of a rhombic diagram T{W) with de, da, and ae tiles. Then 

Fefi(W,ri) F'£fi{W,T2) 

From the above, the definition below is well-defined. 

Definition 2.8. Let LF be a word in {d, a, e}, and let T be an arbitrary tiling of r(fF). Then the 
weight of a word W is 

weight(fF) = ^ wt(F). 

F£{i{W,T) 

^We define the flips more precisely in Section In The well-known fact that one can obtain any tiling of r(TF) from 
any other tiling with a series of flips is proved in |3. 
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In |9], it was shown that the rhombic alternative tableaux satisfy the same recursions as in the 
Matrix Ansatz, and therefore the weight generating function for the tableaux provides the stationary 
probabilities of the process. 

Theorem 2.9 ([^ Theorem 3.1). Let W he a word in {d, a, e}” that represents a state of the two- 
species PASEP with exactly r a’s. Let 

Zn,r = weight (ITQ 

w 

where W ranges over all words in {d, a, e}” with exactly r a’s. Then the stationary probability of 
state W is 

Pr(lT) = —)*—weight(lT). (2.1) 

^n,r 


3 Matrix Ansatz proof of Theorem 2.9 


In this section we give a new proof of Theorem 2.9 by explicitly dehning matrices D, A, and E and 
row vector (n| and column vector \w) that satisfy the hypotheses of a slightly more general Matrix 
Ansatz, and also have a combinatorial interpretation in terms of the RAT. 


3.1 Definition of our matrices 

Our matrices are inhnite and indexed by a pair of non-negative integers in both row and column, 
D {i,j){u,v)]i,j,u,v>o: ^ and E [^{i,j){u,v)\i,j,u,v>0' Our vectors are also 

indexed by a pair of integers, so (n| = [v(ij)]i,j>o and Ira) = [w(u,v)]u,v>o- 

We dehne 'r(ij) = 1 for i = 0,j = 0, and 0 for all other indices. We dehne = 1 for all indices. 

and 0 for all other indices. 


for 0 < u < i and 0 for all other indices. 

U—1 

-I- a[j]q) + 

w=0 

for 0 < u < i and 0 for all other indices. (Here [j]q = -|- ... -I- 1.) 

Since {i,j) specify the row of the matrices, and {u,v) specify the columns, multiplication is dehned 
as 

u,v 

Note that in the case of the matrices H, A, and E, all products are given by hnite sums, since the 
matrix entries are 0 for tt > i -|- 1 or r; > j -|- 1. 







a 
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To facilitate our proof, we provide a more flexible Matrix Ansatz that generalizes Theorem |2 .1 1 with 
the same argument as in an analogous proof for the ordinary PASEP of Corteel and Williams |Tl 
Theorem 5.2], For a word W € {d,a,e}^ with r o’s, we define unnormalized weights f{W) which 
satisfy 

Pr(W) = f{W)IZn,r 

where Zn,r = 'Thw fi^') where the sum is over all words W' of length 

Theorem 3.1. Let X be a constant. Let (tel and |n) be row and column 
D, E, and A be matrices such that for any words X and Y in {d, a, e}, 
satisfied: 

1. {w\XiDE - qED)Y\v) = X{w\XiD + E)Y\v), 

11. {w\X{DA - qAD)Y\v) = X{w\XAY\v), 

III. {w\X{AE - qEA)Y\v) = X{w\XAY\v), 

IV. fi{w\XD\v) = X{w\X\v), 

V. a{w\EY\v) = A(n;|y|n). 

Let W = Wi... Wn with Wi G {d, a, e} for 1 < i < n represent a state of the two-species PASEP of 
length n with r “light” particles. Then 

1 " 

^ {w\A^\v) +^'^{Wi=a) +El^w,=e)\v)- 

' ' ' i=l 


n and with r a’s. 

vectors with (u)||n) = 1. Let 
the following conditions are 


Proof. The proof of Theorem 3.1 follows exactly that of [U Theorem 5.2], Note that the above 
implies that 

_ r r^ HiP + y^ + ET\v) 

□ 


3.2 Combinatorial interpretation of the matrices in terms of tableaux 

Let W be an arbitrary word in {d, a, e} with rhombic diagram r(VP) with the maximal tiling Tmax, 
and let weight(iy) be the weight generating function for Ii{W,Tmax)- Define a free d-strip to be 
a d-strip that does not contain a fi. We call a de or da tile free if the d-strip adjacent to its east 
d-edge is a free d-strip. Note that any de or da tile that is not free must be empty. 

Definition 3.2. Let VP be a word in {d, a, e}. We define the maximal tiling of r(VT) to be the 
tiling that that does not contain an instance of a minimal hexagon (of Figure]^ (b)), for instance 
the tiling of the rhombic diagram in Figure We refer to such a tiling by Tma:^ Tmax can be 
constructed by placing tiles from Pi (VP) inwards, and placing a de tile with first priority whenever 
possible. In other words, all the north strips of Tmax are, from bottom to top, a strip of adjacent 
de boxes followed by a strip of adjacent ae boxes. (A minimal tiling of r(VP) is correspondingly 
defined as the (unique) tiling that does not contain an instance of a maximal hexagon.) 

^Tmax is the unique maximal tiling from [^. 
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Figure 8: Adding a (a) d, (b) a, or (c) e to the end of W. 


We will show that the previously dehned matrices D, A, and E represent the addition of a d- 
edge, an a-edge, and an e-edge to the bottom of r(iy) to form the rhombic diagram T{Wd), 
r(VFa), and T(We) respectively. Recall that these matrices have rows indexed by the pair (i,j) 
and columns indexed by the pair {u,v). We let i represent the number of free d-strips in a tableau 
F G fi(fF, Tmax{W)), and j the number of a’s in W. For the columns, we let u represent the number 
of free d-strips in a tableau F' G fi{Wd,FmaxiWd)) (and respectively, T(yVa) and r(VFe)), and v 
the number of o’s in VFd (and respectively, Wa and We). 

Definition 3.3. Let X{W) be a word in {D,A,E} representing a product involving the matrices 
D, A, and E, corresponding to the 2-PASEP word W in the letters {d, a, e}. 

For example, if VF = deaae, then X(W) = DEAAE. 


Theorem 3.4. Let W be a word in {d, a, e}, and let X = X{W). Then: 


• is the generating function for all ways of adding \W\ new edges of type W to the 
southwest boundary of a rhombic alternative tableau with i free d-strips and j a-strips, to obtain 
a new rhombic alternative tableau with u free d-strips and v a-strips. 

• is the generating function for rhombic alternative tableaux of type W, which have 
u free d-strips and v a-strips. 

• (r(;|A|n) is the generating function for all rhombic alternative tableaux of type W. 


We prove Theorem 3.4 with the following lemma, which says that the matrices D, A, and E are 
“transfer matrices” for building rhombic alternative tableaux with the maximal tiling. 


Lemma 3.5. For the matrices D, A, and E, 


• is the generating function that represents the addition of a d-edge, 

• is the generating function that represents the addition of an a-edge, and 

• is the generating function that represents the addition of an e-edge 

to the southwest corner of a rhombic alternative tableau with the maximal tiling with i free d-strips 
and j a-strips, resulting in a rhombic alternative tableau with the maximal tiling with u free d-strips 
and j a-strips. 


Proof. We describe the possible rhombic alternative tableaux that arise from the addition of a d- 
edge, an o-edge, and an e-edge respectively to the southwest corner of an existing RAT of shape W 
with the maximal tiling, and i free d-strips and j a-strips. 
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The addition of the d-edge to r(iT) does not affect the interior of the tableau, as in the example 
of Figure (a), and the tiling of the new tableau is clearly still a maximal one. Thus for any 
F G ^{W,Tmax{W^))-, we obtain F' G ?i{Wd,Tmax{W^d)) whose weight simply increases by a, the 
weight of the new d-edge. We have thus wt(T') = awt(T). Moreover, the addition of the d-edge 
adds exactly one free d-strip to F. Thus we obtain the desired entry in the matrix D. 

The addition of the a-edge and a vertical strip of adjacent da tiles to the left boundary of 
'TmaxiW) results in a maximal tiling of T{Wa), as in the example of Figure]^ (b). Let us consider 
the entry (i, j), {u,j + 1) of ^ for 0 < u < i. Each free da tile contains either a g or a /3 with no 
restrictions on their positions, for a total of i — u /3’s and u q’s. Thus there are precisely ways 
to choose such a filling of the new tiles. Every such filling contributes a weight of Wa now 

has j + 1 a’s, and it is clear that all other entries of A are zero. Thus we obtain the desired entry 
in the matrix A. 

The addition of the e-edge and a vertical strip of adjacent de tiles followed by j adjacent ae tiles 
to the left boundary of Traax{W) results in a maximal tiling of T{We), as in the example of Figure 
(c). Let us call this strip of new tiles the new e-strip. There are three possible cases for this new 
e-strip. For the following, let us consider the entry (i, j), {u,j) of E for 0 < u < i. 

Case 1: the new e-strip does not contain an a. Then each of the j ae tiles must contain a g, and 
each of the i free de tiles contains either a g or a /3, with no restrictions on their positions, with 
exactly i — u /3’s and u g’s. This gives a total weight contribution of (*)/3*~“g“^L 

Case 2: the new e-strip contains an a in one of the ae tiles. Then each of the ae tiles below that a 
must contain a g, and each of the free i de tiles contains either a g or a /3, with no restrictions on their 
positions, with exactly i — u /3’s and u g’s. This gives a total weight contribution of (^)a/3*~“g“[j](j. 

Case 3: the new e-strip contains an a in one of the free de tiles. Then exactly i — u of the free de 
tiles below the a must contain a /3, and u of them contain a g. This gives a total weight contribution 

Thus we obtain the desired entry in the matrix E. □ 

Proof of \3.4\ The first point is immediate from Lemma 

The second point is due to the following: (tej is a row vector for which the entry with index (0, 0) 
is 1, and the rest are 0. By the first point, ((^c|X)(q o), is, in particular, the generating function 
for adding |LF| new edges of type W to the southwest boundary of a trivial RAT of size 0, to result 
in a RAT of type W with the maximal tiling with u free d-strips and v a-strips. 

The third point is due to the following: |r;) is a column vector with every entry equal to 1. By the 
second point, the generating function for all possible RAT in ^{W,Fmax{W)) is the sum of RAT of 
type W over all choices for the number of a-strips and free d-strips in the fillings. In other words, 
it is the sum over all {u,v) of It follows that {w\X\v) is the desired generating 

function. □ 


3.5 


3.3 Combinatorial proof that our matrices satisfy the Matrix Ansatz 


Using Theorem 3.4, we provide simple combinatorial proofs that our matrices satisfy the equations 


of Theorem 3.1 Let IF be a word in {d, a, e} with r(IF) its rhombic diagram. In this subsection. 
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when we say “addition of a d (or a or e) to W”, we mean adding a d-edge (or a- or e-edge) to the 
southwest point of r(hF), as described in the preceding subsection. 

I. DE - qED = al3{D + E) 

By our construction, consecutive addition of a d and an e to VF results in a de corner with a de 
corner tile as the bottom-most tile of the e-strip that contains it (as well as the right-most tile of 
the d-strip that contains it). This de corner tile contains an a, /?, or q. 

• If the de corner tile contains an a, then the rest of the e-strip containing this tile must be 
empty. Thus the entire e-strip has weight a/3, and the rest of the tableau has the same weight 
as if the de were replaced by a d (with the same hlling in the corresponding tiles). 

• If the de corner tile contains a /3, then the rest of the d-strip containing this tile must be empty. 
Thus the entire d-strip has weight a/3, and the rest of the tableau has the same weight as if 
the de were replaced by an e (with the same hlling in the corresponding tiles). 

• If the de corner tile contains a (/, then this tile has no effect on the rest of the tableau which 
has the same weight as if the de were replaced by an ed (with the same tiling and hlling), and 
the tile itself has weight q. 

Combining the above cases, we obtain that DE = qED + al3{D -|- E), as desired. 

II. DA - qAD = aPA 

By our construction, consecutive addition of a d and an a results in a da corner with a da corner 
tile as the right-most tile of the d-strip that contains it. This de corner tile contains a j3 or q. 

• If the da corner tile contains a /3, then the rest of the d-strip containing this tile must be 
empty. Thus the entire d-strip has weight a/3, and the rest of the tableau has the same weight 
as if the da were replaced by an a (with the same hlling in the corresponding tiles). 

• If the da corner tile contains a q, then this tile has no effect on the rest of the tableau which 
has the same weight as if the da were replaced by an ad (with the same tiling and hlling), and 
the tile itself has weight q. 

Combining the above cases, we obtain that DA = qAD + ajSA, as desired. 

III. AE - qEA = ajdA 

Definition 3.6. We call an ae strip the region of the rhombic diagram that corresponds to a 
maximal a-strip together with an adjacent maximal e-strip. (By maximal a- and e-strips, we mean 
a- and e-strips as they would appear in a maximal tiling of a rhombic diagram, i.e. a strip of adjacent 
da tiles for the a-strip as in Figure]^ (b), and a vertical strip of adjacent de tiles followed by a strip 
of adjacent ae tiles for the e-strip as in Figure]^ (c).) We allow any valid tiling for the ae strip, 
and we call an ae strip maximal if it has the maximal tiling, and we call it minimal if it has the 
minimal tiling. Note that a minimal ae strip has an ae corner tile in the ae corner. 

By our construction, consecutive addition of an a and an e results in a maximal ae strip. For our 
proof, we consider the corresponding minimal ae strip. We apply a series of hips to convert the 
maximal ae strip to a minimal ae strip, and we consider the contents of its ae corner tile. This 
ae corner can contain an a or q. If the ae corner tile contains an a, then the rest of the e-strip 
containing this tile must be empty. Thus the entire e-strip has weight a/3, and the rest of the tableau 
has the same weight as if the e-strip were removed entirely. This operation is the same as if in the 
original tableau, the ae were replaced by an a (with the same hlling in the corresponding tiles). 
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For the other case, if the ae corner tile contains a g, then this tile has no effect on the rest of the 
tableau. Thus the weight of the tableau with the exception of the ae corner tile is the same as the 
weight of a tableau with the same tiling and filling with the ae replaced by an ea. Moreover, this 
new tableau (with the ae corner tile removed from the minimal ae strip) is in fact the maximal 
tableau that corresponds to replacing the ae by an ea. Thus we have as desired, AE = qEA + al3A 
from these two cases. 

Remark 3.7. It is also possible to directly compute the (i, j), {u,v) entry of each term of the 
equations of Theorem |3.H and show that equality holds in each case. 


3.4 Properties of the RAT and enumeration 


Definition 3.8. A flip is an involution that switches between a maximal hexagon and a minimal 
hexagon, and is the particular rotation of tiles that is shown in Figure]^ 

It is a well-known result for rhombic tableaux that one 
can get from any tiling T to any other tiling T' with a 
series of flips. (This is elaborated upon in |9].) In [9l 
Definition 2.16], we extend the flips on a tiling T of r(VF) 
for some word W to weight-preserving flips on the RAT 
T G fi(VF, T), by explicitly defining the weight-preserving 
flip for each possible case of a filling of a hexagon in E. 

This leads to the following definition. 

Definition 3.9. Let LF be a state of the two-species PASEP and T and T' be some tilings of r(LF). 
A RAT E G fi(lT, T) is equivalent to a RAT F' G &{W, T') if F can be obtained from F' by some 
series of weight-preserving flips. 

Let D” be the set of states of the two-species PASEP of size n with exactly r “light” particles. Let 
'I'(n,r) be the set of equivalence classes of RAT whose type belongs to D”. More precisely, ^ G 
is some set of RAT of a single type such that for any F,F' G '(/’) F and F' are equivalent. Moreover, 
if F and F' are equivalent and F € “ip and F' ^ ip', then ip = ip'. 

Finally, from j9], we also have the following theorem. 

Theorem 3.10 (|9] Theorem 2.19). 

n—1 

JJ(a -I- /3 -I- ia/3). 

i=r 


Zn,r{a,l3, 1 ) = 


n 



e e 


Figure 9: A flip from a maximal (left) 
to a minimal hexagon (right). 


This implies the following corollary. 

Corollary 3.11. 

I 'b(n,r) 


/'n\{n + 1)! 
[rj (r + 1)!' 


4 /c-species PASEP 

We now describe a generalization of the two-species PASEP to a fc-species PASEP (also called k- 
PASEP) of a similar flavor. In our new model, we consider k particle species of varying heaviness on 


11 









a one-dimensional lattice of size n. We call the heaviest particle a d particle, followed by ai > 02 > 
• • • > afc-i- For easier notation, we also introduce another particle which we call an e particle to 
represent a hole, and we allow this to be the lightest particle in our set of species. Thus, in our model 
every location on the lattice contains exactly one particle, out of the set of species {d, ai,..., ak-i, e}. 
Moreover, the d particle is allowed to “enter” on the left at location 1 by replacing an e particle at 
that location (with rate a), and it is allowed to “exit” on the right at location n by being replaced 
with an e particle at that location (with rate /3). The particles of type a* are not allowed to enter 
or exit, so we fix the numbers of particles of those species to be r* for i = 1,... ,k — 1. 

For two particle types A and B, we write A > B (respectively, A < B 01 A = B) to mean that A is a 
heavier particle type than B (respectively, A is lighter than B, or they are equal). The dynamics in 
the bulk are the following: a heavier particle of species A can swap with an adjacent lighter particle 
of species B with rate 1 if M is to the left of B, and with rate 0 < qab < 1 if A is to the right of B. 
This means that heavier particles have a tendency to move to the right of the lattice. Our notation 
is shown in the table below: 


A B 

Qab 

d ai 1 < i < k — 1 

qoi 

d e 

^Ooo 

ai aj l<j<i<k — 1 

Qij 

ai e 1 < i < k — 1 

Qioo • 


More precisely, our process is a Markov chain with states represented by words of length n in the 
letters {d,ai,... ,ak-i,e}. The transitions in the Markov chain are the following, with X and Y 
arbitrary words. 

XateY ^ XeaiY XdeY ^ XedY XdaiY ^ XaidY XaiUjY ^ XajaiY 

Qioo QOqo QOi Qij 

eX ^dX Xd^ Xe 


for 1 < i < k — \ and 1 < J < i. 

where hy X ^Y we mean that the transition from X to Y has probability n being the length 
of X (and also Y). 

Definition 4.1. For a given fc-PASEP, we fix n to be the size of the lattice and r* to be the number 
of particles of species Oj for 1 < i < k — 1. We define to be the set of words of length n 

in the letters {d, oi,..., ak-i, e} with r* instances of the letter a* for each i. We also define 

u . 

Remark 4.2. In Section we will provide a Matrix Ansatz solution for the model with different 
parameters qi for every type of transition. However, so far we only have nice combinatorics when all 
the Qi’s are set to equal a single constant q. Furthermore, it is easy to see that it k = 2, we recover 
the 2-species PASEP that we described in the previous section, and if A: = 1, we recover the original 
PASEP. 
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5 Matrix Ansatz solution for the /c-PASEP 


Building on a Matrix Ansatz solution for the usual PASEP by Derrida at. al. |5] and a more general 
solution for the two-species PASEP by Uchiyama in m, we have the following generalization for 
the /c-species process. 

Theorem 5.1. Let W = Wi... Wn with Wi G {d, oi,..., ak-i,e} for 1 < i < n represent a state of 
the k-species PASEP in Suppose there are matrices D, Ai,...,Afc_i, and E and a row 

vector (u)! and a column vector |u) (with (rc||u) = 1) which satisfy the following conditions 


DE — q^ooED — D-\-E, DAi — qQiAiD — Ai, AiE — qiaoEAi — Ai, AiAj — qijAjAi — 0, (5.1) 

(5.2) 


{w\E=-{w\, D\v) = ^\v), 
a p 


then 


^ n k—1 

Prob(Ty) = —-(wl '^(Wi=d) D + l{Wi=e) ^ + X] '^iWi=ai) A 

^n. r-, ri. i . —f 


Tl 


i=l 


where the coefficient of i‘ 


in 


{w\{D + yiAi H-h Vk-iAk-i + EY\v) 

{w\Al^p( ■ ■ ■ A\-\v) 


Proof. For W a word of length n, we define the weight /n(tP) = (^(^| 0^=1 '^(Wi=d) D -|- l(u/i=e) ^ + 

YliZl '^{Wi=aP A\v). 

We show that /n(tP) satisfies the detailed balance conditions 

fn(W) Yi Pr{W^V)= Yi fn{X)Pr{X ^W) (5.3) 

w^v x^w 


for each W G D”', where by Pr(VP — >■ V) and Pr(A — >■ W) we denote the probabilities of the 
transitions W ^ V and X ^ W respectively. This would imply that the stationary probability of 
state W is proportional to fn{W), which would complete the proof. 


We observe that for fixed W, the only terms /^(W) Pr(A — )> V) for some X.,V G D” appearing in 
(7.1), are precisely the terms: 


. fn{eW2...Wn)a, 

• fn{Wi...Wn-ld)P}, 


. and {/,(Wi... Wi-iBCWi +2 ... Wn)A,-fniWi ... W^-iCBWi +2 • • • WnYqBc} where W,W,+i 
BC ioi B > C over 1 < i < n — 1. 


This is because these terms are precisely the terms out of which possible transitions can occur to go 
into or out of W. Moreover, whether those terms appear on the left hand side of Equation (7.1) or 
the right hand side is determined by whether WjlPi+i = BC or IPjVPi+i = CB for B > C. In other 
words, the terms in the bulk are given a sign of (—l)^i'^i+i>Wi) fgj. each i, and the boundary terms 
are given a sign of (—l)^(wi=d) and ( —l)^(Wn=e) fQj- left and right boundaries, respectively. 
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Thus Equation (7.1) can be rewritten as the following: 


^iw,=d or ...Wn) 

+ ^iW„=d or . . . Wn-ld) 

n—1 

+ ^ (/ 4 IT 1 . . . W^.iBiCiWi+2 ...Wn) 

i=l 

- qs^cMWi... Wi.iCiB,W^+2 ... Wn)) 

(5.4) 


where in the above we use Bi := max(Wj, VEi+i) and Cj := min(tTi, ITj+i). 


The reduction rules of Equation (5.1) or ( |5.2[ ) apply whenever VEi = d or e, or ITn = d or e, or 
whenever Wi 7 ^ VEj+i for 1 < z < re. We obtain the following. 


fniW'deW") - qoocfniW'edW") = fn-iiW'dW") + fn-i{W'eW"), (5.5) 

fniW'daiW”) - qoiUW'aidW'') = fn-i{W'aiW''), (5.6) 

UW'aieW'') - qiooUW'eaiW") = fn-i{W'aiW''), (5.7) 

UW'deW”) - qijUW'edW'') = 0, (5.8) 

aU^W”) = fn-i{W"), (5.9) 

lifn{W'd) = fn-l{W'). (5.10) 


For W = Wi ... Wn, we introduce the notation fn_i{W) = /„_i(lTi .. .Wi... Wn) to be the weight 
of the word W with the letter Wi cut out. With this notation, using the reduction rules of Equation 
(5.5), Equation (5.4) becomes the sum oq + oi + ... + an-i + o-n, where 


fLiiW) VTi = 
-fLiiW) VTi = d’ 


/"_i(lE) W„ = d 
-/^i(lE) Wn = e^ 


and 


fn-liW) + fi,t\iW) 



fn-liW) 

-fnt\{W) 

fnt\{W) 

-Pn-liW) 


if WiWi+i = de or ed 
if WiWi+i = dai 
if WiWi+i = ad 
if WiWi+i = aie 
if WiWi+i = eai 


for 1 < z < re — 1 . 


(5.11) 


Notice that for all z > j, the terms fn{W'aiajW") — qijfn{W'ajaiW") = 0. 


Suppose there are a total of s transitions in the bulk. For j = 1,..., s, label the location z where the 
j’th transition occurs (i.e. the j’th z for which Wi d ^i+i) by Wty The strategy of our proof is to 
show that all the fn-i terms that arise from the transitions at the locations cancel with 

other terms Equation (5.11) with an opposite sign. We describe these cancellations in the cases that 
follow. 


(a.) WtjWtj+i = de, so the contribution of terms from this transition is fn-pW) + fn-i (W). Then 
^tj+ibbzj+i+i is necessarily either ed or eat for some t, in which case it contributes the term 
-fnPliW). Similarly, Wt,_,Wt^_ is necessarily either de or o^e for some u, in which case 
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it contributes the term However, the former of these cancels with the term 

and the latter cancels with as desired. 

There are two exceptions to the above. First, if j = 1, then there is no tj-i term. However, 
in this case, W necessarily begins with a d, and so the term cancels with the left 

boundary term —/^_i(VF). Second, if j = n, then there is no tj+i term. However, in this case, 
W necessarily ends with an e, and so the term cancels with the right boundary term 

(b.) WtjWtj+i = ed, so the contribution of terms from this transition is — /*Ti^(bF). 

Then is necessarily either de or ate for some t, in which case it contributes the 

term Similarly, +1 is necessarily either de or da„ for some u, in which 

case it contributes the term However, the former of these cancels with the term 

—fn-i{W), and the latter cancels with as desired. 

There are two exceptions to the above. First, if j = 1, then there is no tj-i term. However, 
in this case, W necessarily begins with an e, and so the term cancels with the left 

boundary term Second, if j = n, then there is no tj+i term. However, in this case, 

W necessarily ends with a d, and so the term cancels with the right boundary term 

/^iW- 

The rest of the cases are similar. Below, we describe the cancellations that occur for each 
transition location. 

(c.) WtjWtj+i = dat, so the contribution of terms from this transition is This term 

cancels with the term since must equal ed or for some u. 

(d.) WtjWtj+i = atd, so the contribution of terms from this transition is This term 

cancels with the term since must equal de or dau for some u. 

(e.) WtjWtj+i = ate, so the contribution of terms from this transition is f^^i{W). This term 
cancels with the term since must equal ed or ea^ for some u. 

(f.) Wt^Wtj+i = eat, so the contribution of terms from this transition is —fn_i{W). This term 
cancels with the term since Wtj_iWt^_^+i must equal de or a^e for some u. 

The cancellations of the boundary terms are treated as the exceptions in cases of (a) and (b). 


It is easy to check from the above that every term cancels with another term in Equatio n (|5.4[ ), 
so indeed, it equals zero. Thus the function fn satisfies the detailed balance in Equation (7.11, as 
desired. □ 


6 /c-rhombic alternative tableaux 

In this section, we introduce a combinatorial object that generalizes the RAT to provide an inter¬ 
pretation for the probabilities of the /c-PASEP. This object, called the A:-rhombic alternative tableau 
(or /c-RAT) is of the same flavor as the RAT, and is similarly defined as follows. 
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6.1 Definition of /c-RAT 


To a word W £ 




we associate a /c-rhombic diagram r(iT) as follows. 


Definition 6.1. Let W G tq be the number of e’s and the number of d’s in 

W. Let an e-edge be a unit edge oriented in the direction —vr. Let a d-edge be a unit edge oriented 
in the direction —7r/2. Let an Oj-edge be a unit edge oriented in the direction —(see Figure 
12). Dehne Pi{W) to be the lattice path composed of the e-, ai-, . .., and d-edges, placed 


end to end in the order the corresponding letters appear in the word W. Dehne P 2 {W) to be the 
path obtained by placing in the following order: tq e-edges, ri ai-edges, r 2 02 edges, and so on, up 
to rk-i Ofc-i-edges, and then d-edges. The fc-rhombic diagram r(LF) is the closed shape that 
is identihed with the region obtained by joining the northwest and southwest endpoints of Pi(W) 


and P 2 {W) (see Figure 10). 


Dehne a lattice path given by W to be composed of the edges in the order they appear in the word 
X, and let us associate this lattice path with the southeast boundary of our rhombic diagram. We 
complete the path to form the diagram by drawing in the following order: to connect the top-most 
corner of the lattice path to its bottom-most corner. 



Figure 10: r(a 2 daiea 2 aieed) dehned by Pi and 
P 2 , with a maximal tiling. 


Figure 11: A 3-RAT of type a 2 daiea 2 aieed of 
weight 


Definition 6.2. A de tile is a rhombus with d and e edges. A dai tile is a rhombus with d and a* 
edges. An Oje tile is a rhombus with a* and e edges. An aiUj tile is a rhombus with a* and aj edges 
for i > j (see Figure [I^. We impose on the tiles the following partial ordering: ajX < aiX' < dX'\ 
and Xe < Xttj < Xui < Xd for i > j and for any edges X,X',X'' . If tile C < tile D according to 
our ordering, we say D is heavier than C. 



Figure 12: e-edge, afc_i-edge, ..., ai-edge, 
d-edge 



Figure 13: A de-tile, doj-tile, Oje-tile, and 
ajOj-tile (with i < j) 


Definition 6.3. A maximal tiling on a Prhombic diagram is one in which tiles are always placed 
from southeast to northwest, and priority is always given to the “heaviest” tiles. 

Dehne a maximal corner to be a corner on Pi(LF) whose edges A and B are such that for any other 
corner on that diagram with edges C and P, AB > CD. The canonical way to tile the rhombic 
diagram with a maximal tiling would be to pick a maximal corner with some edges A and B, and 
place an AB tile adjacent to that corner. The rest of the surface would then itself be a rhombic 
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Figure 14: (a) e-strips and (b) d-strips. 


diagram with the same P 2 . We proceed to tile that surface in the same manner until the untiled 
region has area zero. It is easy to see that such a construction results in a maximal rhombic tiling 
of the fc-rhombic diagram. Let us call this tiling T(W). 


Definition 6.4. An e-strip is a maximal strip of adjacent tiles whose edge of adjacency is an e- 
edge, as in Figure 14 (a). A d-strip is a maximal strip of adjacent tiles whose edge of adjacency is 
a d-edge, as in Figure 14 (b). (This definition is the same for the A:-RAT as it is for the RAT). 

We now define a filling of T{W) with a’s and /3’s as follows. 


Definition 6.5. A filling of a A:-RAT is defined by the following rules. 


• A de-tile is allowed to be empty or contain a or /3. 

• A dai tile is allowed to be empty or contain /3, for each i. 

• An Oje tile is allowed to be empty or contain a, for each i. 

• An ttiUj tile must be empty, for each i > j. 

• Any tile in the same e-strip and above an a must be empty. 

• Any tile in the same d-strip and left of a /3 must be empty. 

Denote the set of fillings of TiW) by fi(VF). We assign weights to a filling F G fi(LF) from the rules 

above by placing a gr in each tile that is not forced to be empty by some a below it in the same 

e-strip, or some (3 to the right in the same d-strip. For an example, see Figure [liP| 

Definition 6.6. Let W G D"', and t be the number of d’s and i the number of e’s in W. For 
F G fi(VF), define the weight wt{F) to be the product of the symbols in the filling of F times 

Define 

^ ^ wt(F) 

W Fefi(W) 

to be the sum of the weights over all A:-RAT corresponding to states in D” ^ Our main result 

for the /c-RAT is the following, which we will prove in the next section. 

Theorem 6.7. Let fi(IT) denote the set of fillings of the rhombic diagram r(VF) with the maximal 
tiling, and let wi{F) denote the weight of a filling in ^(W).Then the stationary probability of state 

®We allow the parameters qsc that represent swapping rates between B-type and C-type particles to vary in 
Section |5| However, to keep the combinatorics “nice”, we fix all these parameters to equal a single constant q. 
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W of the k-PASEP is 


1 




wt(F). 


Conjecture 6.8. Let T' be any tiling of the rhombic diagram r(VL) associated to a state W of 
the /c-PASEP. Let fi(LP, 7*0 denote the set of hllings of tiling T'. Then the stationary probability of 
state W of the fc-PASEP is 

^- Y 


Remark 6.9. When k = 2, the conjecture above is true, as we can dehne a weight-preserving 
bijection between hllings of the arbitrary tiling T' and hllings of the maximal tiling T, in terms of 
“hips” [9j. For k > 3, hips admit the same weight-preserving bijection, but it is no longer necessarily 
the case that we could obtain any tiling T' of a A:-rhombic diagram via hips from the maximal tiling. 


6.2 Matrix Ansatz proof for the fc-RAT 


We will prove Theorem 6.7 using the same strategy as in 


Section for the RAT. 


We provide matrices D, E, Ai ,..., A^-i that correspond to the addition of a d-edge, e-edge, or Oj- 
edge for l<z<A: — Ito the bottom of the path corresponding to a word W of length re to form a 
new rhombic diagram with a maximal tiling of size re -|- 1 that corresponds to the word Wd (or We, 
or Wtti for 1 < i < k — 1 respectively). For A = aft, we show that these matrices satisfy the Matrix 
Ansatz relations 


DE—qED = \{D+E), DAi—qAiD = AAj, AiE—qEAi = AAj, AiAj = qAjAi for i > j. 

( 6 . 1 ) 

The fc-species Matrix Ansatz of Theorem 5.1 would then imply that the steady state probabil¬ 
ity of A:-PASEP state W is proportional to a certain matrix product {w\X{W)\v) with the matrices 
{D, E,Ai,..., Afc_i}. (As in Section]^ we let X(W) be the word in the matrices {D, E,Ai,..., A^-i} 
that corresponds to the word W in the letters {d,e,ai,..., afc_i}.)|^ Similarly to Section 3.2 we show 
that these matrices give a combinatorial interpretation to the construction of the A:-RAT. Therefore, 
the hllings with a’s, /3’s, and g’s of the maximal tilings of the /c-rhombic diagrams provide the steady 
state probabilities for the /c-PASEP. 


In these matrices, the rows are indexed by the tuple {i,ji, ■ ■ ■ ,jk-i) where i is the number of free 
d-strips in a tableau F of the maximal tiling of r(lT) and ji is the number of afs in W. The columns 
of the matrices are indexed by the pair {ifj'i, ■ ■ ■ ,jk-i)i where k is the number of free d-strips in a 
tableau F' of the maximal tiling of r(lTd) (and respectively, TfWe) and r(VPas) for each s) and 
is the number of afs in ITd (and respectively. We and WOg for each s). 


Analogously to the construction of the matrices in the two-species PASEP case, we have now 


and 0 for all other indices. 


D 


(hiivdfc-ihi+ijivjA-i) 


1 

is 


5.1 


in the 


^In Equation ( |6.l[ l, the constant A = ajS is used to slightly generalize the Matrix Ansatz of Theorem 
same manner that Theorem |3.1| generalizes Theorem|2.1[ The statement of the theorem and the proof are very similar 


to that of Theorem 13.11 so we do not provide them here. 
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k-l 

n' 

5=2+1 


for 0 < u < i and 0 for all other indices. 


12—1 

+ a[j]q) + aJ2 

w=0 

for 0 < n < i and 0 for all other indices, where we define j = X^sZiJsi and [j]q = + ... + 1. 

The relations DE — qED = D + E, DAi — qAiD = A*, and AiE — qEAi = + are satisfied by the 
same argnments as in the two-species PASEP case, except with some additional powers of q in the 
eqnations. It remains to show that AtAg = qAgAt for t > s. 

First we compnte the ... Jg,... Jt, ■ ■ ■ ,jk-i){u,ji,... ,jg + 1,... ,jt + I, ■ ■ •, jfc-i) entry of 

AtAg. (The entries of all other indices are antomatically zero). 









a 


i 

W=U 


i 

= E 

w=u 



k-l 

r=t+l 




k-l 

•g + 

r=s+l 




w 


,w-\-u oi—u 


qW+u/^ 


k—1 k—1 

n n (^- 2 ) 


r=t+l 


r=5+l 


Similarly for AgAt, 

i 

— Jl J. + l. is + l,.. Jt + + l jfc-l) 


e(.i)«"' 9‘-" ri n 

r=5+l ^ ^ r=i+l 

k-l k-l 


E n n (6.3) 


r=i+l r=s+l 


It is clear that AfAg = qAgAt, as desired. 


7 A Markov chain on the RAT that projects to the two-species 
PASEP 


We restate here the meaning of a Markov chain that projects to another, and describe the RAT 
as a Markov chain that projects to the two-species PASEP. Snch resnlts exist for the alternative 
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tableaux which project to the regular PASEP, which were originally described in terms of permuta¬ 
tion tableaux (which are in simple bijection with the alternative tableaux) in |2]. Our Markov chain 
has the same flavor as the existing Markov chain defined by Corteel and Williams. The following 
definition is from El Definition 3.20] 

Definition 7.1. Let M and N be Markov chains on finite sets X and T, and let / be a surjective 
map from X to Y. We say that M projects to N if the following properties hold: 


• If xi,X 2 G X with ProbM{xi —> X 2 ) > 0, then ProbM{xi —> X 2 ) = Probiq{f{xi) —>■ /(X 2 )). 


• If yi and 7/2 are in Y and Probi\f{yi —)• 2 / 2 ) > 0, then for each xi G X such that f{xi) = 
yi, there is a unique X 2 £ X such that f{x 2 ) = 2/2 and ProbM{xi —> X 2 ) > 0; moreover, 
ProbM{xi -t X2) = ProbNiyi -t 2/2)- 


Furthermore, we have the following Proposition 7.2 which implies Corollary |7.3| below. 


Let Probm(xo — t x;t) denote the probability that if we start at state xq at time 0, then we are in 
state X at time t. From the following proposition of [2], we obtain that if M projects to N, then a 
walk on the state diagram of M is indistinguishable from a walk on the state diagram of N. 


Proposition 7.2. Suppose that M projects to N. Let xq € X and yo,yi G Y such that /(xq) = 2/o- 
Then 

ProbAr(2/o -t 2/1) = X] ProbM(xo -t xi) 

x' S.t. f(x')=yi 

Corollary 7.3. Suppose M projects to N via the map f. Lety G Y and let X' = {x G X | /(x) = y}. 
Then the steady state probability that N is in state y is equal to the steady state probabilities that M 
is in any of the states x G X'. 


In our case, N is the two-species PASEP (which we call 
the PASEP chain), and M is the Markov chain on the 
RAT (which we call the RAT chain). 


Recall that Qjf denotes the states of the two-species 
PASEP of size n with exactly r “light” particles. We spec¬ 
ify the states of the RAT chain to be the set of 

the RAT equivalence classes of size (n,r), based on the 
fact that different tilings can be chosen to yield equivalent 


tableaux, as mentioned in Remark 3.9 


Now, we define the transitions on T 


{n,r) 


in the RAT chain 


that correspond to transitions on D” in the PASEP chain. 


We introduce the following terminology, as in Figure 15 



Definition 7.4. A corner is a pair of consecutive d and 
e, d and a, or a and e-edges on the boundary of a RAT. 

If there is a de tile, a da tile, or an ae tile (respectively) adjacent to the corresponding edges of the 
boundary, we call that tile a corner tile. 


An inner corner is a pair of consecutive e and d, a and d, or e and a edges on the boundary of a 
RAT. 

®The results in this section could be extended to the fc-RAT in the natural way, but we omit the details and proof 
in this paper. 
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An empty e-strip corresponds to an e-edge on the boundary of the RAT that coincides with its 
top-most boundary. 

An empty d-strip corresponds to a d-edge on the boundary of the RAT that coincides with its 
left-most boundary. 

Lemma 7.5. Let ip G a RAT equivalence class and let F £ ip. If F has a corner of type 

de, da, or ae, then there exists an equivalent F' £ fj that has, respectively, a de tile, a da tile, or an 
ae tile at that corner. 



Figure 16: If the boundary of the tableau contains consecutively d, a, and e, and there is no da tile 
adjacent to the da corner, then a “stack of boxes” as in (a) must occur in the tiling, for some value 
of j. After performing j flips, the configuration in (b) is obtained, with a da tile adjacent to the da 
corner, as desired. 


Proof. First, it is clear that any tiling of a rhombic diagram with a de corner must have a de tile at 
that corner, so for the de case the lemma is obvious. 


Now, for the do and the ae cases, it suffices to prove the lemma for only one of them, since by taking 
the transpose of a tableau and swapping the roles of a and /3, we end up exchanging the d’s with 
the e’s (and consequently the da corners with the ae corners), and so by symmetry, these cases will 
have the same properties. Thus we will prove the da case. 


First, if the da corner already has a da tile adjacent to it, we are done. Thus we assume there is 
not a da tile, which means the tiling of the rhombic diagram must contain the tiles shown in Figure 
16 (a). More precisely, as seen in the figure, the tiles must be a row of j > 1 de tiles on top of j ae 


tiles, with one adjacent da tile on the left. Now it is easy to check that with j flips, we end up with 
the configuration in Figure 16 (b), and moreover, there will a /? in the corner da tile in the tiling (b) 
if and only if there is a /3 in the right-most de tile in the tiling (a) (and otherwise there will be a g). 
Thus with j flips, we obtain an equivalent tableau with a da tile in the da corner, as desired. □ 


Based on the above lemma, we make the following definition: 

Definition 7.6. Let T be a tableau with a corner. We call that corner a q-corner (or an a-corner, 
or a /3-corner) if a tableau T contains a q in the tile adjacent to that corner (or respectively, an a, 
or a /3) for some T that is equivalent to F and has a corner tile adjacent to the corner. 

Definition 7.7. Let T be a tiling of a rhombic diagram F. A d-path on T is a path from some 
point on Pi{F) to some point on P 2 {F) consisting of A- and e-edges. An e-path on 7~ is a path 
from some point on Pi{F) to some point on P 2 {F) consisting of d- and a-edges. We introduce the 
operation of compressing a d-strip in T to obtain a new tiling T' with a d-path in place of the 
d-strip (respectively, e-strip and e-path). We also introduce the inverse operation of blowing up 
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Figure 17: Let di and d 2 be the indicated d- and e-paths on T- Then T is the compression of 
bl(T, di) and bl(T, ^ 2 ) at the highlighted d- and e-strips, respectively. 


a d-path in T' to obtain a new tiling T" with a d-strip in place of the d-path (respectively, e-path 
and e-strip). 

Compressing a d-strip means selecting its northern border to be the d-path, and then gluing together 
the north and south e-edges and a-edges of each tile in the d-strip, thereby replacing the d-strip by 
the d-path. Similarly, compressing an e-strip means selecting its western border to be the e-path, 
and then gluing together the west and east d-edges and a-edges of each tile in the e-strip, thereby 
replacing the e-strip with the e-path. If s is a d- or e-strip of T, then we denote by com(T, s) the 
new tiling T' that results from compressing at s. 


For the inverse, blowing up a d-path means replacing each e-edge of the path with a de tile, and 
each a-edge with a da tile, to obtain a d-strip from the new tiles. Similarly, blowing up an e-path 
means replacing each d-edge of the path with a de tile, and each a-edge with an ae tile, to obtain 
an e-strip from the new tiles. If p is a d- or e-path of T, then we denote by bl(T, p) the new tiling 
T” that results from blowing up at p. Figure 17 illustrates these dehnitions. 


By convention, if p is a path of length 0, then blowing up p results in replacing it by an empty 
e-strip or an empty d-strip (depending on whether p coincides with the west boundary or the north 
boundary of the rhombic diagram, respectively). Conversely, compression of an empty e-strip or an 
empty d-strip results in replacing those strips with a single point. 


It is easy to see that compressing is the inverse of blowing up. 

Let T be a RAT of size {n,r,k) with tiling 7~, and let 'ipiF') G ^(n,r) denote the equivalence class 
that F belongs to. Below we describe the RAT chain transitions on F, which are also transitions 
on 'ijj{F). 


7.1 “Heavy” particle enters from the left. 

If F has an empty e-strip e, then there is a transition in the RAT chain from F that corresponds to 
a “heavy” particle entering from the left in the PASEP. Let the type of F be eW. 

We dehne a new RAT T as follows. Let p be the south-most point on Pi{F) (the southeast boundary 
of F) such that there are exactly n — k — 1 e- and a- edges on Pi{F) southwest of p. Let p be any 
d-path originating at p. Let T' = bl(com(7~, e), p). It is easy to check that T' is a valid tiling of 
r(dX) which has size (n, r,k + 1). 

If n — — 1 > 0, the new d-strip of T' is non-empty, so we place a /3 in its right-most tile, which 

is valid since that tile must be either a de tile or a da tile. Furthermore, p was chosen to be the 
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(b) Wd We 


Figure 18: For both examples, let the left tableau have tiling T, and denote the indicated empty 
(e- or d-) strip by e and the marked (e- or d-) path by p. Then we obtain a new diagram with 
tiling bl(com(T, c),p), and (a) a /3 is placed in the resulting d-strip or (b) an a is placed in the 
resulting e-strip. 


south-most point such that there are n — k — 1 e- and a-edges southwest of it, so the right-most 
tile of the new d-strip is also the bottom-most tile of the a- or e-strip it lies in, and thus does not 
interfere with the rest of the filling of the tableau. We define ProbRAr(.F —?■ T) = The weight 

of F with the exception of e equals the weight of T with the exception of the newly added d-strip. 
The weight of the new d-strip of T is afi, and the weight of e is (3. Therefore, wt(T) = ^wt(T), 

and so wt(F) Frob rat{F T) = 

For the exceptional case, if n — A; — 1 = 0, then the newly added d-strip of T is empty, and thus 
has total weight a. In this case, the PASEP state corresponding to F is of the form ed^~^, and the 
PASEP state corresponding to T is d*^. Then wt(F’) = wt(T) = and so in this case 

we have wt(F) ProbRAT(E —)• T) = . 

7.2 “Heavy” particle exits from the right. 

If F has an empty d-strip e, then there is a transition in the RAT chain from F that corresponds to 
a “heavy” particle exiting from the right in the PASEP. Let the type of T be VFd. 

We define a new RAT T as follows. Let p be the east-most point on Pi{F) such that there are 
exactly r -|- A: — 1 d- and a- edges on Pi{F) northeast of p. Let p be any e-path originating at p. Let 
F' = bl(com(T, e), p)- It is easy to check that F' is a valid tiling of T{We) which has size (re, r. A: — 1). 

If r -|- A: — 1 >0, the new e-strip of F' is non-empty, so we place an a in its bottom-most tile, which 
is valid since that tile must be either a de tile or an ae tile. Furthermore, p was chosen to be the 
east-most point such that there are r -|- A: — 1 d- and a-edges northeast of it, so the bottom-most tile 
of the new e-strip is also the right-most tile of the a- or d-strip it lies in, and thus does not interfere 
with the rest of the filling of the tableau. We define ProbreAr(E —)• T) = The weight of F 

with the exception of e equals the weight of T with the exception of the newly added d-strip. The 
weight of the new d-strip of T is a/3, and the weight of e is a. Therefore, wt(T) = — wt(T), and so 

wt(F) ProbRAr(E ^ T) = 
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For the exceptional case, if r + /c — 1 = 0, then the newly added e-strip of T is empty, and thus 
has total weight /3. In this case, the PASEP state corresponding to F is of the form and the 

PASEP state corresponding to T is e”. Then wt{F) = wt(T) = and so in this case we 

have wt(F) Proh rat{F T) = 




(c) da —^ ad for a /3-corner 



Figure 19: (a) and (b) show transitions at an a-corner, and (c) and (d) show transitions at a 
/3-corner. 


7.3 A “heavy” particle exchanges with a “hole”. 

If F has a de corner, then there is a transition in the RAT chain from F that corresponds to a 
“heavy” particle swapping places with a “hole” in the PASEP. Let the type of F be WdeY, and 
suppose it has tiling T. The de corner necessarily corresponds to a de tile. This tile contains an a, 
a /3, or a q. We describe these three cases below. 


7.3.1 The de corner tile contains a /3. 


We define a new RAT T as follows. Let the d-strip containing the de corner tile have length A. Let 
p be the south-most point on Pi{F) such that there are exactly A — 1 e- and a- edges on Pi{F) 
southwest of p. Let p be any d-path originating at p. Let F' = bl(com(T, c),p). It is easy to check 
that T' is a valid tiling of r(lTedy), as in Figure 19 (d). 


If A — 1 > 0, then we place a /3 in the right-most box of the newly inserted d-strip 5. Such a 
filling is valid since the right-most box (containing the new /3) is necessarily the bottom-most box 
of the e- (or a-) strip that contains it, and so s does not interfere with any of the other tiles in 
T. We define Probi{Ar(E —)• T) = The weight of T equals the weight of F. Therefore, 

wt(F) ProbijAr(E ^ T) = 
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If A — 1 = 0, then necessarily F corresponds to a PASEP state WdedP for some j, and T corresponds 
to the state WeF~^^. The newly added d-strip is empty, and so wt{F) = /3wt(T). Therefore, 
wt(F) PiohRATiF ^T) = 

7.3.2 The de corner tile contains an a. 

We define a new RAT T as follows. Let the e-strip containing the de corner tile have length A. 
Let p be the east-most point on Pi{F) such that there are exactly A — 1 d- and a- edges on Pi{F) 
northeast of p. Let p be any e-path originating at p. Let T’ = bl(com(T, e),p)- It is easy to check 


that T' is a valid tiling of r(lTedy), as in Figure 19 (b). 


If A — 1 > 0, then we place an a. in the bottom-most box of the newly inserted e-strip s. Such 
a filling is valid since the bottom-most box (containing the new a) is necessarily the right-most 
box of the d- (or a-) strip that contains it, and so 5 does not interfere with any of the other tiles 
in T. We define Probi{Ar(T —)• T) = The weight of T equals the weight of F. Therefore, 

wt(F) Prob,jAr(T ^ T) = 


If A — 1 = 0, then necessarily F corresponds to a PASEP state e^deY for some j, and T corresponds 
to the state e'^"*“^dy. The newly added d-strip is empty, and so wt{F) = awt(T). Therefore, 
wt(F) PiohRATiF ^T) = 

7.3.3 The de corner tile contains a q. 

We define a new RAT T by simply removing the de corner tile from F. We define ProbijAT(T —?• T) = 
Since a single tile of weight q was removed, wt(T) = (7 wt(r). Therefore, wt(T) ProbijAr(T 
_ gwt(r) 

^ ) ~ JV+1 ■ 

7.4 A “heavy” particle exchanges with a “light” particle. 


If F has a da corner, then there is a transition in the RAT chain from F that corresponds to a 
“heavy” particle swapping places with a “light” particle in the PASEP. Let the type of F be WdaY. 


By Lemma 7.5, we can assume that F has a da tile at the da corner. This tile contains a /3 or a g. 


We describe these two cases below. 


7.4.1 The da corner tile contains a /3. 


We perform exactly the same operation as for the de case containing a /3. Once again, we define 
PiohRATiF —)• T) = ]y^. In all but the exceptional case, the weight of T equals the weight of F. 

_ wt(T) 

N+1 ■ 


V+l ■ 

Therefore, wt(F) PiohRATiF 


T) = 


In the special case, if F corresponds to a PASEP state Wdad^ for some j, and T corresponds to the 
state ITad-^"''^, then we have wt(F) = /3wt(r). Therefore, wt(F) ProbijAT(T —t T) = . 
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7.4.2 The da corner tile contains a q. 


We perform exactly the same operation as for the de case containing a q. Again, wt{F) ProbRyi'7’(T —>• 

rr\ — gwt(T) 

^ ) ~ N+l ■ 


7.5 A “light” particle exchanges with a “hole”. 


If F has an ae corner, then there is a transition in the RAT chain from F that corresponds to a 
“light” particle swapping places with a “hole” in the PASEP. Let the type of F be WaeY. By Lemma 


7 . 5 , we can assume that F has an ae tile at the da corner. This tile contains an a or a q. We describe 
these two cases below. 


7.5.1 The ae corner tile contains an a. 

We perform exactly the same operation as for the de case containing an a. Once again, we define 
Probij^7’(T —)• T) = In all but the exceptional case, the weight of T equals the weight of F. 

Therefore, wt{F) Pioh rat{F T) = 

In the special case, if F corresponds to a PASEP state e^aeY for some j, and T corresponds to the 
state then we have wt(F) = awt(T). Therefore, wt(T) ProbRAT(-E —)■ T) = . 

7.5.2 The ae corner tile contains a q. 

We perform exactly the same operation as for the de case containing a q. Again, wt(T) ProbRAr(E —>■ 

rji\ _ (j'Wt(r) 

^ > ~ Af+l ■ 

7.6 A lighter particle type exchanges with a heavier particle type. 

We describe only the WedY —>■ WdeY transition, but the same holds true for WadY —>■ WdaY and 
WeaY —7- WaeY if the corresponding letters are used. If F has an inner ed corner, then there is 
a transition in the RAT chain from F that corresponds to a “hole” swapping places with a “heavy” 
particle in the PASEP. Let the type of F be WedY. Then to form the tableau T, we simply append 
a de tile to the outside of F, adjacent to the ed inner corner. We place a q inside the tile, and thus 
obtain a valid filling T of type WdeY with a (7 in its de corner. 

We define ProhjiATiF —^ T) = Therefore, since (7wt(F) = wt(T), we have wt(F) PiohRAriF —> 

'T^\ _ wt(X’) 

^ > ~ “/v+T* 

The operator pr is clearly a surjective map from the set to It is easy to see by our 

description of the transitions on the RAT chain that it indeed projects to the PASEP chain. 

7.7 Stationary probabilities of the RAT chain 

We carefully summarize the transitions out of a RAT F (and consequently from the equivalence 
class of F), depending on the chosen corner at which the transition occurs. We will be referring 
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Figure 20: Some of the transitions on some of the states in All the transitions involving the 
circled tableaux are included. 


to these cases further on. First we make the following definitions. Let F have size (n, r, k) and let 
A = (Ai,..., \k+r) be the partition given by the lengths of the d-strips from top to bottom. Assume 
that A has at least one non-zero part. 


Definition 7.8. We define Xr be the indicator that equals 1 if F has an empty e-strip, and 0 
otherwise. We define Xl be the indicator that equals 1 if F has n empty d-strip, and 0 otherwise. 


Definition 7.9. We call a q-corner a corner that contains a q. (Refer to Definition 7.6 for the 
precise definition.) We call a top-most corner an a- or /3-corner such that the length of the d-strip 
containing it equals Ai. (If the corner in the top-most position contains a q, we do not call it a 
top-most corner). We define the indicator 6^ which equals 1 if the top-most corner contains a /?, 
and 0 if it contains an a. Analogously, we call a bottom-most corner an a- or /3-corner such that the 
length of the row containing it equals the length of the smallest non-zero row of A. (If the corner 
in the bottom-most position contains a g, we do not call it a bottom-most corner). We define the 
indicator which equals 1 if the bottom-most corner contains an a, and 0 if it contains a /3. We 
call a middle eorner an a- or /3-corner that is neither a top-most corner or a bottom-most corner 
(and not a ( 7 -corner). 


7.7.1 Summary of transitions F — )• T 

Denote by 7 r(F —)• T) the rate of transition from tableau F to T (where by rate we mean the 
unnormalized probability). We obtain the following cases for the transitions from F to T. 

1. For a transition at a middle corner, a top-most corner with (3^ = 1, or a bottom-most corner 
with 5^ = 1, we have wt(r) = wt(F), and tt{F —)• F) = 1. 

2. For a transition at a top-most corner with 6^ = 0 such that the length of the e-strip containing 
it is greater than 1, we have wt(T) = wt(F) and tt{F —>■ T) = 1. Then the top-most corner of 
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T will be an a-corner. 

3. For a transition at a bottom-most corner with (5^ = 0 such that the length of the row containing 
it is greater than 1, we have wt(T) = wt(-F) and 7r(-F —)■ T) = 1. Then the bottom-most corner 
of T will be a /3-corner. 

4. For a transition at a top-most corner with 6^ = 0 such that the length of the e-strip containing 
it is 1, we have wt(T) = ^ wt{F) and tt{F —)■ T) = 1. 

5. For a transition at a bottom-most corner with = 0 such that the length of the d-strip 
containing it is 1, we have wt(T) = ^ wt{F) and tt{F —)■ T) = 1. 

6 . For a transition at an empty e-strip, we have wt(T) = Q!wt(F’) and Tr{F T) = a. T will 

not have an empty e-strip, and it will have a top-most corner that contains a (3. 

7. For a transition at an empty d-strip, we have wt(T) = f3wt{F) and Tr{F T) = /3. T will 

not have an empty d-strip, and it will have a bottom-most corner that contains an a. 

8 . For a transition at an inner corner, we have wt(T) = qwt{F) and Tr{F T) = q. 

9. For a transition at a g-corner, we have wt(T) = ^ wt(T) and tt{F —T) = 1. 

Our main theorem is the following. 

Theorem 7.10. Consider the RAT ehain on il' the RAT equivalenee elasses of size {n,r). Fix 
a RAT F and its equivalenee elass V’- Then the steady state probability of state tp is proportional to 
wt(F). 

Proof. To prove the theorem, it suffices to show that for each RAT F, the following detailed balance 
condition holds. Let TZ be the set of RAT such that there exists a transition from F to T G TZ. Let 
S be the set of equivalenee elasses of RAT such that for each G S, there exists some S G such 
that there is a transition from S to F. Though we actually work with the equivalence classes, we 
write for simplicity S G S. 


wt(F) = wt{S)7r{S ^ F). (7.1) 

Ten ses 

Let the RAT F have type W. First we treat the transitions going out of T to T G TZ. By 
the construction of the RAT chain, it is clear that there is a transition with probability 1 for 
every corner (including the top-most-, bottom-most-, middle-, and g-corners), a transition with 
probability a for an empty e-strip, a transition with probability /3 for an empty d-strip, and a 
transition with probability q for every inner corner. These transitions directly correspond to all the 
possible transitions out of the two-species PASEP state W. Suppose F has Cq (/-corners, C a- or 
/3-corners, and I inner corners. Thus we obtain 

^ ^ Tr{F -G T) = C -\- Cq -\- ql -\- crdi + fddn. C^-^) 

TGR 

For the transitions going into F from some S G S, we observe that any transition from one tableau 
to another ends with a (/-corner or an a- or /3-corner, an empty e-strip, an empty d-strip, or an inner 
corner. Thus it is sufficient to examine all such properties of F to enumerate all the possibilities for 
S G S. We examine the pre-image of the cases for the possible transitions going into F to obtain 
the following cases for S. 
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5. 


1. For a middle corner, a top-most corner with 6^ = 0, or a bottom-most corner with 6^ = 0, we 
have wt(S') = wt(F’) and 7r(S' F) = 1. This is the inverse of Casej^of Section 7.7.1 This 
gives a contribution of wt(T)(C' — 2 -|- (1 — 6^) -|- (1 — (5^)) to the right hand side (RHS) of the 
detailed balance equation]^ 

2. For a top-most corner with 5^ = 1 and 6ji = 0, we have a transition involving an empty e-strip 
of S, so wt(5) = ^ wt(F) and Tr{S F) = a. This is the inverse of Case of Section 


7.7.1 


This gives a contribution of wt{F)6^{l — Sr) to the RHS of the detailed balance equation. 

For a bottom-most corner with 6^ = 1 and Jj;, = 0, we have a transition involving an empty 
d-strip of S, so wt(5) = ^ wt(F) and tt{S F) = /3. This is the inverse of Casej^of Section 
This gives a contribution of /34 wt(T)<5^(l — 6l) to the RHS of the detailed balance 


7.7.1 


ion. 


equa 

For a top-most corner with 6^=1 and (5/^ = 1, there are two possibilities. For the first, S 
could fall into Case of Section 7.7.1 meaning that the top-most corner of S' is a /3-corner, 
which results in the usual transition with wt(S) = wt(F). For the second possibility, S could 
fall into Case 1^ of Section 7.7.1 meaning that the top-most corner of S is an a-corner and 
the column containing it has length 1. In that case, wt(S) = awt{F). In both situations, 
7 r(S —?• T) = 1. We obtain a contribution of wt{F)5j^ {5 r + a{l — 5 r)) to the RHS of the 
detailed balance equation. 

For a bottom-most corner with 6^ = 1 and 6l = 1, there are two possibilities. For the first, 
S could fall into Case of Section 7.7.1 meaning that the bottom-most corner of S is an 
a-corner, which is the usual transition with wt(S') = wt{F). For the second possibility, S 
could fall into Case of Section 7.7.1 meaning that S has a bottom-most corner containing a 
/3 and the row containing it has length 1. In that case, wt(5) = f3wt{F). In both situations, 
7r{S —>■ F) = 1. We obtain a contribution of wt{F)6^ {6 l -|-/3(1 — Sr)) to the RHS of the 
detailed balance equation. 


For a g-corner, we have wt(5) = ^ wt(F) and 7 r(S' F) = q. This is the inverse of Case of 
Section 7.7.1 We obtain a contribution of wt(F) to the RHS of the detailed balance equation. 


F) = 1. This is the inverse of Case 
We obtain a contribution of gwt(F) to the RHS of the detailed balance 


For an inner corner, we have wt(S') = gwt(F) and 7 r(S' 
1^ of Section 7.7.1 


equation. 


We sum up the contributions to the RHS of the detailed balance equation to obtain 


wt(S)7r(S 

SeS 


F) = wt(F)(C + Co + qI-5f-5!^ + 6^(1 - 6 r ) + d^(l - 6 l ) 


R/ 


+ 5f{5R + «(1 - 5r)) + 5^{5l + /3(1 - 5l))). (7.3) 


We see that after simplification. Equation 7.3 equals Equation |7.2[ so indeed the desired Equation 


7.1 holds for “most” F, save for the easily-verified degenerate cases. 


□ 


The proof above circumvents the use of the Matrix Ansatz, and is another way to prove our main 
result of Theorem 12.91 

®Note that if C < 2, the formulas we give have some degeneracies. However, it is easy to verify that these do not 
cause any problems due to cancellation of all the degenerate terms. 
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As an example, we show some of the transitions on the RAT chain for tableaux of size (4,1) in 
Figure!^ 
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